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Abstract 

We propose a scheme for the realization of a quantum walker and a quantum simulator for the 
Dirac equation with ultracold spinor atoms in driven optical lattices. A precise control of the 
dynamics of the atomic matter wave can be realized using time-dependent external forces. If the 
force depends on the spin state of the atoms, the dynamics will entangle the inner and outer degrees 
of freedom which offers unique opportunities for quantum information and quantum simulation. 
Here, we introduce a method to realize a quantum walker based on the state-dependent transport 
of spinor atoms and a coherent driving of the internal state. In the limit of weak driving the 
dynamics is equivalent to that of a Dirac particle in 1+1 dimensions. Thus it becomes possible to 
simulate relativistic effects such as Zitterbewegung and Klein tunneling. 

PACS numbers: 03.75.Mn,03.67.Ac,03.65.Pm 
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I. INTRODUCTION 



Quantum simulators aim at the simulation of complex quantum systems in well control- 
lable laboratory experiments Such a simulation is especially useful when the original 
quantum system is experimentally not accessible and classical simulations are impossible 
due to the exponential size of the Hilbert space. Furthermore, quantum simulators offer the 
possibility to tune the experimental parameters to explore novel physical phenomena. Im- 
portant examples include the simulation of solid state systems with ultracold atoms j^J, Dirac 
dynamics with graphene or trapped ions flfl and sonic b.ack ho.es in Bose-Einstein 
condensates [7j. 

Ultracold atoms in optical lattices are especially suited for such a task, since their dy- 
namics can be controlled with an astonishing precision and their dynamics can be measured 
in situ. In the present paper we propose to use spinor atoms in tilted or driven optical lat- 
tices to realize a quantum walker - a paradigmatic system in quantum information science 



12[ - and a quantum simulator for relativistic Dirac dynamics. In contrast to previous 



proposals which were based on the realization of artificial gauge fields using atoms with a 
tripod internal structure [3, Q , we focus on simple spinor atoms with two internal states, 
which are used routinely in ongoing experiments. The necessary correlations between the 
internal and the external dynamics can be realized with state dependent external forces, in 
particular by magnetic gradient fields 
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and microwave transitions between the internal states 
the full dynamics of the atoms. 
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e driving of this external field 



21] . one can obtain control over 



II. BLOCH OSCILLATIONS AND TRANSPORT IN OPTICAL LATTICES 

A Bloch oscillation is the counter-intuitive dynamics of a quantum particle in a periodic 
potential subject to a static external force. The force accelerates the particle until it reaches 
the edge of the Brillouin zone, where it is Bragg reflected. If Landau-Zener tunneling to 
higher bands can be neglected, this leads to a fully periodic motion. This peculiar kind o: 
dynamics was postulated by Bloch already in 1928 in the context of electrons in crystals '± 
but never observed because of the strong scattering of the electrons. Thus it took another 
seven decades before Bloch oscillations could be demonstrated for electrons in semiconductor 
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FIG. 1: (Color online) Dynamics of ultracold Cesium atoms in a driven optical lattice: (a) A static 
external field leads to Bloch oscillations, (b) A sinusoidal driving of the field leads to directed 
transport. We assumed an optical lattice with a wave length A = 1064 nm, a depth of Vq = En 
and a maximum field strength of Fq/tti = 0.42 m/s 2 . The time is given in units of the Bloch period 
Tb = 8.44 ms and the period of the external driving Td r i ve = 10 ms, respectively. 



superlattices [231 ] . photons in waveguide arrays [24J and ultracold atoms in optical lattices 



251 ] . Especially the latter realization shows an astonishing level of coherence and offers 



various possibilities to precisely control the atomic dynamics. Recent experiments demon- 
strated Bloch oscillations over tenth of thousands of periods and macroscopic distances \}\ . 
An example of Bloch oscillations of an atomic matter wave is shown in Fig. [T] (a). 



Extensive possibilities to control the atomic motion can be rea 



ized in driven optical lat- 



tices, i.e. in lattices with a time-dependent external field |15l.l26N31|. A particular interesting 
case is a periodic driving, when the direction of the field is reversed before the matter wave 
reaches the edge of the Brillouin zone, decelerating the wave packet back to a standstill. 
Thus the atoms always have a positive momentum such that directed transport is realized, 
as illustrated in Fig. [T] (b) for a sinusoidal driving. 

To be precise, we have simulated the dynamics of ultracold Cs atoms in an optical lattice 
with a wavelength of A = 1064 nm and a depth of Vq = E^. The strength of the external 
field was assumed to be Fn/ra = 0.42 m/s 2 , which can be easily realized by accelerating the 

n nn 

optical lattice |25[ or by a magnetic gradient field [15|, [17j. The dynamics of the atomic 



matter wave is then given by the Wannier-Stark Hamiltonian 

H{t) = + V cos 2 (k x) + F{t)x, (1) 

where fco = 27r/A and Er = f^k^jlm is the atomic recoil energy. The initial state has been 
chosen as a Bloch state with zero momentum weighted by a Gaussian with a width a = 6 A. 
For a static field one then finds the celebrated Bloch oscillations with a period given by 
the Bloch time Tg = AttH/XFq. The maximum displacement of the wave packet is given by 
d = A/Fq, A being the width of the ground Bloch band 

A directed transport of the atoms can be realized by means of a time-periodic driving. 
During one period of the driving the atoms are first accelerated and the decelerated back to a 
standstill, such that the atomic wavefunction is displaced in space but otherwise unaffected. 
An example for directed transport in a driven optical lattice is shown in Fig. [1] (b) for a 
sinusoidal driving 

F(t) = F lC os(27rf/T), (2) 

assuming a driving strength of Fi/m = 0.42 m/s 2 with a period of T = 10 ms. The transport 
properties of a driven optical lattice can be calculated analytically within the tight-binding 
approximation 
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27l |. It is found that transport is possible only for a pure ac-driving, 



while it is generally forbidden in the case of a combined dc- and ac-field F(t) except for the 
case of a resonant driving T = uTb- The actual transport velocity depends on the initial 
state, but an upper bound can easily be found: 



la \ noj 

where J n is an ordinary Bessel function. The case of a pure ac-driving corresponds to n = 0. 
The dispersion of a wave packet is generally negligible for an initially broad wave packet as 



it vanishes with l/a 4 la, |27|, |28| . This approach to controlling transport in driven optical 



lattices was experimentally demonstrated in 30l. 131]. 



III. QUANTUM WALKS 

Now consider the dynamics of a spinor atom, where the atoms experience a different 
field strength depending on their initial state. Using Bloch oscillations and the directed 
transport described above, it is possible to entangle the position and the internal degrees 
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FIG. 2: (Color online) Quantum walk of ultracold Cs atoms in a driven magnetic gradient field, 
(a) The coin operations: Rabi frequency of the microwave 7r/2-pulses inducing transitions between 
the two hyperfine levels, (b) The shift operation: A driven magnetic gradient field induces state- 
dependent transport, (c) The resulting dynamics of the total atomic density \tp^(x, t)\ 2 + \ip±(x, t)\ 2 . 
Parameters are given in the text. 

of freedom. This provides the basis for the implementation of a quantum walk or a Dirac 
quantum simulator. 

One possibility to induce a state dependent transport of spinor atoms is the use of a 
magnetic gradient field. In particular we consider the dynamics of ultracold cesium atoms 



with the internal states | 7) = \F = A,m F = 4) and | \) = \F = 3,m F = 3) 16]. 
The magnetic moments of these two internal states are opposite, such that they move into 
different directions in a magnetic gradient field. The effective potential is given by 

V mF (x) = g F m F fi B AB x, (4) 

assuming a linear magnetic gradient field B z (x) = ABx. The effective Lande factor of the 
two internal states is given by g F =4 = 1/4 and g F =z = —1/4 32| and fis = 9.274 x 10~ 24 J/T 
denotes the Bohr magneton. The magnetic gradient field can also be varied in time to control 



the atomic motion 



15] . Furthermore, we assume that the internal state of the atoms can be 



manipulated coherently by resonant microwave pulses with Rabi frequency Q(t) 20|. The 
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FIG. 3: (Color online) Characterization of a Bloch quantum walk: (a) Increase of the standard 
deviation Ax 



(x) 2 . (b) The total atomic total atomic population \ip-f(x, t)\ + \ipi(x, t)\ e 



after 15 driving periods. The parameters are given in the text. 

dynamics of the atomic state (i/)^(x,t),il)±(x,t)) is thus given by the Hamiltonian 

-h 2 d 2 



H{t) 



+ Vq cos (kox) 



+ 



m(t) 



(5) 



2m dx 2 

V t (x) 
V^x) 

A coherent quantum walker can be implemented by a sinusoidal driving of the gradient 
field AB(t) = ABq x sin(27rt/T). During one period, the atom are first accelerated and then 
decelerated back to a standstill as illustrated in Fig. [1] (b). However, the transport direction 
depends on the direction of the external field, which is opposite for both internal states, 
thus realizing an effective quantum walk. The coin operation is then realized by a ir/2 pulse 
coupling the two hyperfine levels after integer multiples of the driving time T. An example of 
a quantum walk is shown in Fig. [2] (c), where we have plotted the dynamics of the complete 
atomic density \ip^(x,t)\ 2 + \ipi(x,t)\ 2 . The time- dependence of the magnetic gradient field 
and the microwave pulses is illustrated in the upper panels of the figure. Initially, all atoms 
are assumed to be in the internal state | \) <g> ipo{ x ), where the spatial wavefunction i/jo(x) is 
a Bloch state with k = 0, weighted with a Gaussian envelope of width a = 3 A. The depth 
of the optical lattice is given by Vq = Er. 

It has to be noted that the magnitude of the force is different for the two spin states leading 
to a trivial overall displacement of the atomic density pattern. Furthermore this introduces 
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FIG. 4: (Color online) Zitterbewegung of spinor atoms in a driven magnetic gradient field, (a) 
Evolution of the position expectation value for different values of the rotation angle 6. (b,c) Evolu- 
tion of the probability densities \ip^(x,t)\ 2 and \ipi(x,t)\ 2 for 9 = O.Itt. The remaining parameters 
are given in the text. 

a phase shift between the two spinor components, which has more severe consequences. This 
phase shift can be avoided by choosing internal states with the same mj?, i.e | t) = |4, 3) 
instead of |4, 4), or it can be compensated actively by a static magnetic field. 



IV. DIRAC DYNAMICS 



For weak driving, the d yna mics of a quantum walker is equivalent to that of a Dirac 
particle in 1+1 dimension |33|], which can be seen as follows. During one period T, the 
atoms are first displaced depending on their internal state and then the internal state is 
rotated by an angle 9. Assuming that the internal rotation is fast compared to T, the 
dynamics is given by the evolution operator U(nT) = (c^fA)™ with 

I Dd \ 

U 1 = and U 2 = exp(-i6a x /2), (6) 

\ D_ d j 

where denotes the translation operator over a distance d. In the limit of weak driving, 
i.e small values of d and 8, the dynamics can be described by an effective Hamiltonian, 
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FIG. 5: (Color online) Dynamics of spinor atoms in the Dirac regime with an additional potential 
step, (a-c) Evolution of the atomic density (left) and the position expectation value (x(t)) (right) 
for different heights of the potential step: V s /Er = (a), 0.015 (b) and 0.1 (c). (d) Position 
expectation value (x(t)) after a fixed propagation time ifi na i = 150 ms as a function of the potential 
step V s . Results of a full quantum simulation (solid blue line) are compared to the effective 
discrete Dirac approximation (open circles) . The remaining parameters are the same as in Fig. 0] 
with 9 = 0.1 ir. 

U(nT) « exp(— iH e sn), which is given by 

6 

H cS = dpa z + -a x . (7) 

This is just the Dirac Hamiltonian with an effective mass m = 9/ 2d and the momentum 
operator p = —id/dx. 

According to these results, ultracold spinor atoms show an effective relativistic dynamics 
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if the displacement d and the rotation angle 9 are reduced. This is simply realized by 
increasing the lattice depth Vo and decreasing the strength of the microwave driving field. 
An example is shown in Fig. H]for atoms initially prepared in the internal state | \) = |3, 3) 
with a width a = 10 A. The optical lattice is now deeper as in the previous examples, 
V = 5E R . Transitions to the other spin state | f) = 1 4, 3) are induced at integer multiples of 
the period Td r i ve with a variable rotation angle 9 = f2i pu i se . For 9 = 0, the internal state does 
not change such that the atoms are steadily transported by the time-periodic external force. 
Increasing 9 and thus the effective mass of the Dirac particle leads to a Zitterbewegung of the 
atoms, which is illustrated in Fig. H](a). As expected, the period of the oscillations decreases 
with the effective mass m. Microscopically, the oscillations of the position expectation value 
(x(t)) result from the interference of the two components of the atomic state, which are 
shown in Fig. H] (b,c) 

One of the most surprising predictions of Dirac theory is the Klein paradoxon: Relativistic 
particles are not repelled by a strong repulsive potential, but perfectly transmitted if the 



height of the potential step exceeds 2mc 2 34J. Signatures of this effect can be readily 



explored with ultracold atoms in driven optical lattices. We consider a situation where an 
additional blue detuned tophat laser beam is focused aside of the atom cloud. We model 
the induced optical dipole potential by a tanh profile such that the total optical potential is 
given by 

V(x) = ^ att ice(z) + V stcp (x) (8) 
= Vo cosmos) + -^(tanh((x — x s )/w) + 1). 

If the additional potential is weak enough, the directed transport mechanism remains mostly 
unaffected. However, the atomic wave packets accumulate a dynamical phase depending on 
their position which is equivalent to a potential in the effective Dirac hamiltonian ([7|). In 
this regime, the dynamics is described by the discrete evolution operator 

U{nT) « (U 3 U 2 Ui) n with 

U 3 (x) = exp(-iV step (x)T/h) (9) 

and Ui t 2 given in Eqn. ([6]). The dynamics in the presence of an additional potential is shown 
in Fig. [5j comparing a full simulation to the Dirac approximation (Q. A good agreement is 
observed for weak potential steps, V s < 0.02 E R . 




FIG. 6: (Color online) Quantum walk in the presence of decoherence. (a) Increase of the standard 
deviation Ax(t) for different values of the dephasing rate k. (b) The standard deviation in a loglog 
plot showing a different scaling Ax(t) ~ t a in the quantum regime (n = 0) and in the classical 
regime (k = 100s _1 ). The dashed black line is a linear fit with the result a = 0.94 and a = 0.59, 
respectively, (c) The standard deviation Ax(tfi na i) after a fixed time ifj na i = 150 ms as a function 
of K. 

A pronounced signature of Klein tunneling observed in Fig. [5] is that the repulsive po- 
tential step enhances the transmission of the atomic matter wave. Without the potential, 
the atoms exibit a Zitterbewegung as discussed above. If a weak additional potential step 
is included, e.g. = 0.015 Er, the oscillating motion is gone and the atoms move steadily 
above the barrier. An even stronger potential, however, suppresses the directed transport 
underlying the Dirac quantum simulator. The effective description by the Dirac hamilto- 
nian ([7]) is no longer appropriate and the atoms are reflected just like ordinary Schrodinger 
particles. In this regime the Dirac approximation clearly deviates from the exact simulation 
results. 
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FIG. 7: (Color online) Atomic Zitterbewegung in the presence of decoherence. Shown is the position 
expectation value (x(t)) for different values of the dephasing rate k. Parameters are the same as 
in Fig. H with 9 = 0.2 vr. 

V. DECOHERENCE 

Real experiments with spinor atoms are limited by a loss of coherence due a coupling to 
the environment. The schemes discussed above make use on the opposite magnetic moment 
of the two hyperfine states | f) and | \) to transport the atoms in different directions. 
However, this makes the atomic state vulnerable to dephasing caused by magnetic field 
fluctuations 9], while spatial coherence can be conserved over seconds To explore the 
effects of decoherence, we simulate the dynamics in the presence of pure dephasing described 
by the Master equation 

f = ~i[H, P\~\ {pip + pal - 2a z pa z ) . (10) 



For the numerical calculations we use the quantum jump method [35J, averaging over 100 
trajectories. 

Decoherence of the atomic states essentially turns a quantum random walk into a classical 
one, which is confirmed by the simulations shown in Fig. [61 After a short initial period the 
standard deviation Ax grows much slower in the the presence of dephasing. The linear 
increase Ax(t) ~ t gradually changes into the classical diffusion law Ax(t) ~ t 1//2 if k 
is increased. The different scaling is analyzed quantitatively in part (b), where Ax(t) is 
plotted on a loglog scale together with a linear fit, omitting the initial expansion of the 
atomic matter wave. The different slope of the two curves directly reveals the different 
diffusion exponents. The fit yields an exponent of a = 0.94 and a = 0.59 for k = and 
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k = 100 s _1 , respectively, which is in good agreement with the theoretical expectation a = 1 
for a quantum and a = 1/2 for a classical random walk. A quantitative analysis of the 
decoherence process is provided in Fig. [6] (c), where the standard deviation Ax after a fixed 
propagation time tg na i = 150 ms is plotted as function of k. One observes that Ax decreases 
with k until the random walk becomes completely classical at k ~ 80s -1 . 

Also the Dirac quantum simulator discussed in the previous section is vulnerable to 
decoherence of the atomic hyperfine states. Figure [7] shows an example of the atomic Zit- 
terbewegung for the same parameters as in Fig. H] for different values of the dephasing rate 
K. The oscillations of the mean position of the atoms are much less pronounced for k > 0. 
We find that the dephasing rate must not exceed 20 s _1 to observe a clear signature of a 
Zitterbewegung. For stronger dephasing no interference effects and thus no oscillations are 
visible. 

VI. CONCLUSION 

In the present paper we have explored the rich dynamics of spinor atoms in tilted and 
driven optical lattices. These systems are nowadays routinely realized experimentally and 
can be controlled with a high accuracy. We have discussed how a state-dependent transport 
can be realized by a periodic driving and how this can be used to implement a quantum 
simulator. A quantum random walk is realized when transitions between the internal states 
are driven by microwave or Raman pulses. In the continuum limit of weak driving, the 
dynamics is given by the Dirac equation in 1+1 dimension. Thus it is possible to investigate 
relativistic effects such as Zitterbewegung or Klein tunneling in a well controllable laboratory 
experiment. 

All elements needed for the implementation of the proposed quantum simulator have been 
demonstrated quite recently such that an experimental realization is possible with current 
technology. A major obstacle can be decoherence of the internal atomic state. On the 
other hand, decoherence processes offer the possibility to experimentally study the quantum 
classical transition for a Dirac particle. 
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